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A transitional layer matching the asymptotically flat exterior of a charged, mas- 
sive toroidal body to an interior with spatially cylindrical symmetry is described. 
The changes in the geometry, which by themselves would require an energy tensor 
violating the energy conditions of classical general relativity, are compensated for 
by an additional strong electric field. Part of its gradient is consumed, depending 
on how much the exterior toroid departs from local cylindrical symmetry; what 
is left over can be used to effect further transitions in the cylindrical interior. 
An example is given of such a transition creating an angular deficit, allowing an 
otherwise infinite cosmic string to be captured within a bounded system. 



1 Introduction 

There exist many exact solutions for cylindrically symmetric spacetimes in 
general relativity. Some of them model physically mundane sources, such 
as linear electric currents. [HE] Others, such as cosmic strings, are more ex- 
otic. [3l HI [5] Regardless of what values their metric components approach 
as r — > co, the infinite extent of their sources in the ±z direction make 
them inadequate as approximations to asymptotically flat spacetimes with 
reasonably bounded sources. Yet the much greater complexity of coordi- 
nate systems adapted to spatially compact bodies which locally approach 
the symmetry of linelikc sources make it difficult to model more realistic 
analogs. 

In many cases, it seems intuitively obvious that the physical properties 
of a successfully modelled system with cylindrical symmetry should carry 
over to the local vicinity of long but finite rods or very thin torii. [6] In 



other cases, where a cylindrical model has some bizarre pathology, those 
global differences are made into an argument against such generalization. 

m 

In the absence of such solutions, it is still possible to adapt cylindri- 
cally symmetric spacetimes as models of bounded bodies by embedding 
them inside a transitional layer that on the outside is shaped like a finite 
body, in this case a thin torus. Because the transitional metric components 
represent real physical changes in the local geometry, this layer will have 
a non- vanishing energy tensor. By themselves, such fields would in places 
violate the energy conditions of classical general relativity. However, by 
including an additional strong field which by itself satisfies the energy con- 
ditions but which does not interfere with the transition from one topology 
to another, the violating contributions of the transitional functions can be 
compensated for so that T ab as a whole satisfies these conditions. 




Figure 1: Cross-section of transitional ring. Toroidal exterior, coordinate 
< r\ < ?7i increases inward from zero at spatial infinity; cylindrical inte- 
rior, coordinate < r < r\ increases outward. External boundary at 770, 
internal boundary at 771 = — ln(ri/>t). Angular coordinate 9 is common to 
both coordinate patches: 9 = at the outer equator and ±ir at the inner 
equator of the ring. The center of the ring is to the left of the figure. 
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The exterior spacetime will be that of a thin, nearly uniformly electri- 
cally charged massive ring of radius k in the background coordinates [8], 
modelled by the 'conformastationary' line element [HI El HTj 

2 2 

ds 2 = a- 2 (dt - Wi dx l f - - ^(dri 2 + d6 2 + sinh 2 ? 7 dip 2 ) (1) 

(cosh?? — cos#)^ 



a(r),6) = l + V(ri,9) (2a) 

V(r],9) = — ^v/cosh rj — cos 9 P_ 1/ / 2 (cosh77) (2b) 
2y 2 

where x % can be any of the spatial coordinates and P_i/ 2 (x) is the Legendre 
polynomial of order —1/2. Both the ADM mass and the absolute value of 
the net charge of the ring, as seen by a distant observer, are nq. Among 
the many advantages of conformastatic spacetimes is their isotropy. A 
consequence of this, unlike the more general Weyl spacetimes, is that their 
physical interpretation is unambiguous: what looks like a ring to a distant 
observer, still looks like a ring to a local observer. 

If Wi is constant it will not appear in any component of T ab , so for 
convenience it will be set to in the remainder of this paper. 

The signature will be [H ] and geometric units used such that cq, 

the speed of light in the limit V — > 0, and the gravitational and electric 
force constants G and k e are all equal to 1. r = (x 2 + y 2 ) 1 ! 2 denotes the 
cylindrical radial coordinate in the interior space. Where r appears as a 
subscript it denotes a partial derivative, f r = df /dr; all other subscripts 
are labels only. 

The definition of standard toroidal coordinates is given in an appendix. 



2 Transitional region 

2.1 Metric 

Consider the line element 

ds 2 T = a~ 2 dt 2 -a 2 K 2 {N 2 e- 2 ^[dTi 2 + d9 2 } + M 2 dip 2 ) (3) 
with, for now, a = 1. N(tj,6) and M(r],6) are piecewise defined by 



N(v,0) = { 



2 sin ( —x 



(cosh i] — cos 9) 

(cosh r] — cos ( 
2 



< 7] < T]o 



rji < rj < oo 
(4) 
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M( v ,e) 



sm — x 

.2 



sinh 77 
(cosh 77 — cos ( 

sinh ri 



where 



(cosh r) — cos 9) 
1 

_ V-Vo 

X = A 

Arj 

A?/ = 771 - r/o 



< 77 < 770 
% < < r?i 

7/! < 7] < CO 

(5) 



For < ?y < 770 ([21) is the Minkowski metric in standard toroidal coor- 
dinates (see Appendix). 77 parameterizes a series of nested torii centered 
on a ring of coordinate radius re, where 77 — > 00; and — ir < 8 < tt and 
< (p < 2tt are the angular coordinates of the minor and major orbits 
respectively. 

For 77 > 771 the coordinate transformations 



r = 

z = 



(6a) 
(6b) 
(6c) 



show the interior spacetime to be also flat, but locally perfectly cylindrical, 
with < z < 2ttk a linear but periodic coordinate, and metric 



ds 2 



2 dt 2 -a 2 (dr 2 + dz 2 



de 2 ) 



(7) 



Examine now the transitional region 770 < r) < Tji, and generalize to the 
case where a(x l ) is a function of the spatial coordinates other than (p. In 
the frame diagonalizing the pure-electric part of T ab (not the same as the 
frame diagonalizing T ab as a whole) 

87ra 4 T aS = diag[l,-l,l,l] | Va| 2 + diag[l, 0, 0, 0] [-2aV 2 a -aVa-VlnM] 

+ derivatives of N, M 

It can be seen that T ab can be split into three parts. The first is 
a pure electric field, which by itself satisfies the energy conditions. Its 
gradient can be adjusted by adjusting the net charge q in (j2b)) and the 
outer boundary 770. The third part is due only to the transition functions. 
Its magnitude decreases as 770 — > 00 as the exterior space more closely 
approaches cylindrical symmetry, for then the relative ^-dependence of M 
and N decreases and M — > 1, N — > 2. 

The middle part contains an a-only part and an interaction term. The 
latter 's relative contribution will also decrease with increasing 770. The 
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Laplacian of a is partly dependent upon the choice of q and ij , but it 
may also be to some extent adjusted independently. By giving a a small 
but adjustable negative second derivative, T 00 can be made positive and 
dominant over any other component of T ab . 

The two primary constraints on a are that it smoothly (C 1 ) matcrQ 
with its counterpart on the exterior and that its gradient at r\\ be positive 
(or zero, if a flat interior is desired). The first constraint may be satisfied 

by 



a{rj,B) = l + (l-nx 2 )(C+Ai-i Dx) I sin 2 -x 



cosh 77 — cos( 



1/2 



COS —X 

2 , 



cosh 77 

where the constants C and D are the values at 770 of 

— —= \J cosh 77 P-1/2 (cosh 77) 
2y 2 

and its first derivative respectively. For 770 > 3 they approach [12j 

C ps (7/ +ln4) g /27r (8a) 
D ps q/2ir (8b) 

The second constraint is satisfied by restricting 

_ DA V _ A77 
MS/W- 2C + 3DAr] ~ 2r7 + 3Ar7 + 2.77 l J 

2.2 Energy conditions 

For static systems, the pointwisc weak, strong, null, and dominant energy 
conditions are satisfied if for the fully diagonalized energy tensor |T3] 

T 66 > (10a) 
T 66 ± jHi > o ^ 10b ) 

1 " ■ \ I " ■ I) (10c) 

Of these seven conditions, the one typically coming closest to violation is 
T° 5 - at 77 = 770, 9 = it. Satisfying this requires that 

7T 2 / 37T 



A* > Mmin W 1 + e ""° (H) 

2 V OTo / 

With the upper limit on /1 fixed by ([5]), which cannot be greater than 1/3 
in any case, (jTTJ) sets the minimum charge nq a ring of outer boundary 770 
must have to effect the transition to the cylindrical interior @ 

1 Propcrly 'smooth' requires the metric be C°° continuous across all boundaries, in 
other words that it match to all orders of its derivatives. 'Smooth' in this paper refers 
only to C 1 continuity, or continuity in the metric itself and its first derivatives. 

2 In preparing the table in section [3] fJ, m i n was calculated numerically rather than 
using the approximate formula JTTJ. 
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Larger values of [i may be chosen to reduce the gradient of a at the 
endpoint rji to zero if desired, so as to match the toroidal exterior to a 
flat and empty interior (section [4. 2p . However, the advantage of keeping it 
as large as possible is that the trick here can be repeated with additional 
transitional layers to alter other metric functions. Not only does the cylin- 
drical symmetry of the interior simplify the math needed to do this, it also 
allows the introduction of features difficult to model in vacuum spacctimcs 
outside bounded sources, such as an angular deficit (section [4. 3p . 

2.3 Optical effects. A gravitational hall of mirrors 

The interior space is truly cylindrical. Geodesies propagating parallel to 
the central axis remain parallel, even though the toroidal exterior has made 
the interior periodic in the z direction. A person looking in the ±z direction 
would sec the back of her own head. At any given r, the speed of light in 
any direction is the same whether one happens to be closer to the 'inside' 
or 'outside' edge of the outer ring, since the interior metric includes no 
dependence on 9. Yet the spacetime as a whole is asymptotically flat. 

3 Limitations 

At 971 , relabeled r\ in the new coordinates, a has the form 



a{ri)\ Vl = l+E + F^-^) (12a) 

E = (1 - fi)(C + D Ar)) (12b) 

F = D(Ar]) 2 ~ fi(Ar)){2C + 3D A V ) (12c) 

which after the transformation (|6aj) becomes 

a(r)\ ri = i+E-F(r-n)/ri (13a) 

= l + E-F\n(r/ri) (13b) 

E = (g/27r)(l-/i)ln(8«;/ri) (13c) 



F = (q/2n) [(1 - /i)(A^) 2 - 2/x(Ar/) ln(8«/ri)] (13d) 

If the gradient of a has not been brought to zero at this point, it may 
be left as (|13b[) . representing the electric field generated by a long wire 
centered on r = with mass and charge per unit length F (section |4. 1 1) . 
Alternatively, additional cylindrical shells transitioning to other interior 
geometries may be placed inside. The maximum possible distortion of the 
internal geometry that may be further effected is primarily determined by 
the relative field gradient at t\ 

C = -(ra r /a)\ ri =F/(1+E) (14) 
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Choosing the minimal attenuation of Va from its external value that 
still satisfies the energy conditions (/_t = /i m i„) and choosing A77 = 1, 
the threshold values of E and Q for selected 770 and q, 5 < 770 < 30 and 
—4 < log 10 q < 4 are given in the table below. 



770 


E/q 


q = 


10- 


4 


10- 2 


1 


10 2 


10 4 


30 


5.15 


1.59 


• 10 


-5 


1.51 • 10~ 3 


2.59 


•HT 2 


3.08- 


10~ 


2 


3.09- 10~ 2 


25 


4.36 


1.59 


• 10 


-5 


1.53- 10~ 3 


2.97 


• lO" 2 


3.64- 


10- 


2 


3.65- 10~ 2 


20 


3.56 


1.59 


• 10 


-5 


1.54- lO" 3 


3.49 


• lO" 2 


4.46- 


10- 


2 


4.47- 10" 2 


15 


2.77* 


1.09 


■ 10 


-5 


1.54- 10~ 3 


4.22 


• lO" 2 


5.73- 


10- 


2 


5.75- 10" 2 


10 


1.97* 








7.88 • lO" 4 


5.30 


• lO" 2 


7.99- 


10- 


2 


8.03- 10" 2 


6.7 


1.44* 










4.89 


• lO" 2 


9.69- 


10- 


2 


9.77- 10~ 2 


5 


1.14 














6.35 • 


10- 


2 


6.51 • 10~ 2 


For the lowest-*; cases 


iii 


this table, E/q = 


2.74, 1.94, 1.43 


for 


1)0 


= 15,10,6.7 



respectively. 

Since \i m i n C 1, £ for each 770 is nearly a constant multiple of g, as 
would be expected from (|13c|) . As 770 decreases and the exterior toroidal sur- 
face departs from local cylindrical symmetry to a greater extent, a greater 
attenuation in Va is required to make the transition and a greater mini- 
mal charge q is necessary to keep fi m in < Mmaxi which is why there are no 
solutions for the lowest 770, q pairs. 

As the charge is increased, £ approaches asymptotic maxima for each 
770 (figure [2]). This maximum increases up to about 770 = 6.7: the field 
gradient at 770 increases as the torus becomes a more compact body; but 
beyond this point, the increasing departure from local cylindrical symmetry 
saps more of the outer field strength in order to effect the transition and 
C, the field strength remaining at 771, again declines. For the choice of 
transition functions used here the maximum possible £ is less than 0.10, 
although minor tweaking of the transition functions can get it above that 
threshold]^ 

4 Examples 
4.1 Charged wire 

Leaving a(r) in the logarithmic form (|13b[) makes the interior spacetime 
vacuum except for a cylindrically symmetric electric field, representing (af- 
ter matching to a suitable second interior solution at some r B < r\) the 

3 For instance, globally replacing sin 2 (^-a;) — > sin 4 ( ^a: ), cos 2 (2-x) — » 1 — sin 4 (^x), 
which reduces the curvature of N and M at rjo where l|10|l come closest to violation. 
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0.08 



0.06 



0.04 — 



0.02 




log(10)q 



Figure 2: Maximum values of ( = —ra r /a at r\ as a function of log 10 q for 
selected 770. 



exterior to a long wire centered at r = with a linear mass and charge 
density Q = F. As noted above, the periodicity of z means that the wire 
is not really infinitely long, although the metric, energy tensor, and most 
local measurements will be identical with those of a wire that is. 



4.2 Flat space 

Truncating Vet — > at some r A < r\ will match the spacetime to a flat 
interior (figure A suitable function matching to a{r\) (eq. (|13ap ) and 
a(r) = a A = constant at r A is 

a(r) = a A {l-vy 2 ) (15) 



2n - (Ar 
V = (r-r A )/(Ar) 
Ar = fx — T A 



(16) 



The interior space is flat and isotropic, though still periodic in the z- 
direction. Spatial distances are greater and time runs slower relative to the 
outside by a factor of a A = (1 — v)^ 1 (1 + E). The local speed of light in 
all directions is (l/a A ) 2 . 
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Flat Flattening Toroidal — > cylindrical 




Figure 3: a(r), a(rj) for the example of section 14.21 A second transi- 
tional shell covering r A < r < r\ is placed within the original toroidal- 
to-cylindrical transitional shell, eliminating the remaining gradient of a(r) 
to create a flat space inside. 



4.3 A captive string 

By transitioning to an interior with an angular deficit, a cosmic string 
with its exact vacuum solution may in effect be captured within a spatially 
bounded system (figure I3J. Write the line element as 

ds 2 = a- 2 dt 2 -a 2 (dr 2 +dz 2 +b(r) 2 d6 2 ) (17) 

with a, Ar, and y defined as in section (|4.2|) and b(r) defined afl 

r > ri 

r A < r < n (18) 
r < r A 

4 Properly since r\ marks the outer boundary of the cylindrical interior space, the 
coordinate r should not be continued beyond it as in figure [4] and equation 118t . This 
coordinate abuse is done for clarity of exposition. 



b(r) 



r + a(l-y) 2 
B(r + r s ) 
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Figure 4: b(r) = y/—goe for the example of section l4~3l The slope of \/—gee, 
db(r)/dr, is attenuated from 1 to a lesser value over r A < r < r\ to create 
an angular deficit within the interior space. Since r' = r + r s , the true 
origin r' = or r = —r s , where the cosmic string is located, is shifted 
inward from r = 0. 



1 - B 
2 

1-B\ 
B 



Ar (19a) 

(19b) 



For B = 1 — 4A < 1 and the coordinate substitution r' = r + r s , the 
spacetime inside r < r A is the exterior to a straight, static cosmic string 
with linear mass density A, together with a radial electric field. 

The general energy tensor has again an electric part T£ b = (Stto 4 )" 1 | Va| 2 
x diag[l, — 1, 1, 1] plus 

T 6fl_ T 66 = (8 7 ra 3 6)" 1 (-26a rr - ab rr - 2b r a r ) (20a) 
T^-Tf = (87ra 3 fe)- 1 a6 rr (20b) 

With a, b, cr, b r , b rr > and <x r , a rr < 0, the energy conditions will be 
satisfied if T 65 - T zi > 0, or 

— a rr > b~ 1 (ab rr + b r a r ) (21) 

Again, the energy deficit will be met by attenuating Va from its value at r\. 
For simplicity, Va will be reduced to at r A as in section 14.21 eliminating 
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the electric field and leaving only the vacuum exterior to a string located 
at r' = 0. In other words, the entire amount of C will be spent on creating 
the largest possible angular deficit. Then (f2"Tj) will be satisfied so long as 
a = ar A v, a < 1. Equating this er to (|19a[) . and choosing a = 1 and 
Ar <C r\ so that r A w r% , gives for the maximum possible angular deficit 



allowing at best a relative reduction of 10% in ^Jgee under the constraints 
of this model. 

Returning to the original form of the exterior metric ((T|) , it can be seen 
that a spinning string, so long as it is rotating around its longitudinal axis, 
can be modelled as well: the introduction of a constant-coefficient term 
we dO has no effect on the energy tensor anywhere, including the transition 
layer and the exterior toroidal spacctime. 

A rotating loop of string could also be modelled with the introduction 
of a constant-coefficient w z dz metric term in the internal coordinates, with 
a corresponding constant-coefficient w v dip in the exterior metric. But then 
the outer spacetime would have a conical singularity along the axis of the 
ring and so would no longer be asymptotically flat. 

In general, the upper limit on the relative field gradient £ is the chief 
restriction on what cylindrical interiors may be embedded within these 
mediating rings. A transitional layer designed to contain a van Stockum 
cylinder |14j . for instance, would require using up part of £ on bringing 
w(r) = —gte/gtt to a non-constant value and another part on creating an 
angular deficit. This could be done for a causally well behaved interior 
solution; but the maximum possible £ < 0.10 is too weak to bring these 
parameters to values that would introduce closed timelike curves. 
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2n( 
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Appendix. Standard toroidal coordinates 

The definitions of standard toroidal coordinates are given here, rather than 
in the main body of the text, to avoid confusion between the cylindrical 
coordinate system (p, Z, tp) in the external space by which they are most 
easily defined and the cylindrical coordinate system (r, z, 9) used through- 
out this paper for the cylindrical space within the toroidal shell. 

With reference to a cylindrical coordinate system in the exterior space- 
time (figure [S]), toroidal coordinates can be implicitly defined by 

P = u a ( 23 ) 

cosh T] — COS 9 

„ k sin 9 

Z = (24) 

cosh r\ — cos 9 

where k is the radius of the source ring. They have a flat three-metric 

ds 2 = K 2 (coshry-cos6>)- 2 [(drj 2 + d9 2 )+ sinh 2 r] dtp 2 ] (25) 

Axisymmetric solutions to Laplace's equation suitable for the exterior to a 
source have the form 

f n (r],6) = (cosh?y — cos0) 1 / 2 P Il _i/2(cosh?7) (AcosnO + BsinnO) 

(26) 

They are built from Legendre polynomials of half-integer order, sometimes 
called 'toroidal harmonics'. These do not have a closed form. For n > 
they become proportional to e~( n+l /*) , i for r\ — ► oo. [H] 

Note that the solutions to Laplace's equation arc no longer simply sep- 
arable (they belong to the 'R-scparable' class of solutions [5]). Even the 
zero-order solution will have some relative dependence on the angular co- 
ordinate 9. But like the (cosh -q — cos 9) term in the metric itself, this 
dependence will rapidly decrease as rj gets large. 

For further discussion of toroidal coordinates, toroidal harmonics, and 
modelling the exterior fields of ring sources, see [51 112). 
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Zk 




► 

p 



Figure 5: Standard toroidal coordinates, constant-y? cross-section. Cen- 
tral ring located at Zr = 0, p = k (black dots). Radial- type coordinate -q 
(constant-?7 surfaces are black lines) increases from on the vertical axis 
and at spatial infinity to oo on the core of the ring; angular coordinate 9 
(grey lines) increases from on the Z = plane outside the ring, to ±7r 
on the central disk spanning it. Adapted from [16] . 
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